Muon diffusion in diluted paramagnetic systems 
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Muon Spin Rotation/Relaxation (/iSR) is a powerful probe to study diffusion of hydrogen- like 
particles. One of the most common model to analyze /iSR signals from materials where diffusion 
occurs, is by using functions derived from the stochastic theory of Kubo and Toyabe along with the 
strong collision model. Unfortunately the current formalism can not be used to analyze data from 
systems where not only muons diffuse, but also the magnetic moments of the material are dynamic. 
In this article we introduce a new model which accounts for both dynamical effects, and use this 
model to analyze /iSR signals from Dyo.4Yi.eTi207. 

PACS numbers: 76.75. +i, 66.30.j- 



I. INTRODUCTION 

The study of muon diffusion is a topic of interest not 
only for science but also for technology, since the muon 
can be thought of as a light proton 1 - 2 . Muon diffusion 
is measured using the Muon Spin Rotation/Relaxation 
(/iSR) technique. In this technique polarized muons are 
implanted in a sample, and the positrons arising from 
their asymmetric decay are recorded as a function of time 
and direction^. With this information one may recon- 
struct the time dependent polarization signal of the muon 
ensemble. The analysis of this signal provides important 
information about the local magnetic environment and, 
if present, about the muon diffusion rate. 

The most common model used to analyze ^SR sig- 
nals from systems where muons diffuse is based in the 
probabilistic theory of R. Kubo and T. Toyabe 4 together 
with the strong collision approximation 4 ^^. In its sim- 
plest version this model assumes that, after the muons 
are implanted, they sense a local static field which is 
random and has a Gaussian distribution. For para- 
magnetic systems this local magnetic field is typically 
originated by the nuclear magnetic moments, which are 
static within the /zSR time window. The model as- 
sumes also that the muons can hop between minima 
of potential with a hopping rate Uh, and therefore the 
muons see instantaneous field changes with this same 
rate. This model, which is referred as the Dynamical 
Gaussian Kubo- Toyabe (DGKT) model, is suited to ana- 
lyze signals from materials where the magnetic moments 
are densely packed. For magnetically diluted systems 
though, the internal field distribution is not Gaussian 
but Lorentzian 7 ^. 

No standard model exist to analyze signals if both, 
muon diffusion and fluctuations of the magnetic moments 
of the system, are present. As we will see in the next sec- 
tion, in magnetically dense systems a DGKT-likc model 
cannot disentangle muon hopping rates from fluctuations 
of the magnetic moments of the system. Nevertheless, 
this is not the case for magnetically diluted materials, 



and we report here a model suited to analyze signals from 
this type of systems. Using numerical simulations of this 
model we have found that magnetically diluted systems 
are specially well suited to study muon diffusion, since 
its onset produces a change on the shape of the relaxing 
signal (even if the hopping rate is much smaller than the 
fluctuation rate of the magnetic moments). After this, 
in Section HU we present a set of phenomenological for- 
mulas derived from our simulated polarization functions, 
and which can be used to analyze experimental data. In 
Section IIVI we use these formulas to analyze data from 
Dyo.4Yi.6Ti 2 07; and then in section |Vl we discuss the 
results and comment on the limits of the model. 



II. A MODEL FOR MUON DIFFUSION WITH 
MAGNETIC MOMENT FLUCTUATIONS 

In the Gaussian Kubo- Toyabe (GKT) models the lo- 
cal field at each muon site is taken randomly from a Gaus- 
sian field distribution: 



P G (B, A) = P G (B X , A) py(B y , A) P^(B Z , A) (1) 



where: 
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cxp(- 7 2 J B 2 /2A 2 ) 



(2) 



In this last equation 7 is the muon gyromagnetic ratio 
(2tt x 135.54MHz/T), A is the width of the field dis- 
tribution, and Bi is one of the three components of the 
magnetic field. Since the local field at each muon site is 
static, then the polarization of the muon ensemble can 
be analytically calculated and is given by: 



A(t) 



(1- A 2 t 2 )exp(-A 2 t 2 /2 



(3) 
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The GKT model describes systems of Heisenberg mag- 
netic moments which fill completely the sites of the lat- 
tice and which are static. If the field sensed by the muons 
changes due to fluctuations of these magnetic moments, 
then the strong collision model can be used to model the 
effect of these fluctuations on the muon polarization. In 
this strong collision models the field change is a Marko- 
vian processes with a characteristic fluctuation rate Vf 
(i.e. after a time the local field changes instanta- 

neously to other taken randomly from the field distribu- 
tion in Eq. The resulting dynamical model is the 
DGKT one mentioned in the introduction. Now, if the 
local magnetic moments are static but the muons are dif- 
fusing at a rate Vh (and assuming that the field at the 
new site where the muon hops is uncorrelated with the 
previous one), then the field sensed by the muons will 
be effectively fluctuating with a rate v^, and the strong 
collision model can be also applied to describe the muon 
diffusion process. Then, in the DGKT model fluctuations 
of the magnetic moments are completely equivalent to a 
muon diffusion process; therefore, in magnetically dense 
systems it is not possible to distinguish one dynamic pro- 
cess from the other. This is not the case for magnetically 
diluted systems as we will show in the following. 

Lets assume a static scenario again, but now most of 
the magnetic moments have been removed (i.e. substi- 
tuted) from the lattice. In this diluted system the muons 
which thcrmalize near to a magnetic moment can sense 
higher fields than those which land far from a magnetic 
moment. This means that many magnetically inequiva- 
lent muon sites are present in a diluted system. Mathe- 
matically, this can be model by introducing a distribution 
of A. This distribution is given by-i£: 



P(A,a) 
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2 exp(-a 2 /2A 2 ) 



(4) 



where a is the characteristic width of the distribution. 
If the field distribution in Eq. [T] is averaged with the A 
distribution on Eq. |4j the total magnetic field distribu- 
tion seen by the muon ensemble is the three dimensional 
generalization of the Lorentzian distribution as required 
for a diluted magnetic system^: 



Pl (B, a) 



7 3 a 



7r 2 (a 2 + 7 2 B 2 ) 2 



(5) 



where B is the magnitude of the local field. The static 
polarization signal which arises from this field distribu- 
tion is: 



1 2 

A(t) = - + -(!- at) exp(-at) 



(6) 



If the field at the muon site changes due to fluctua- 
tions of the magnetic moments of the material, then the 
dynamical behavior is introduced with the strong colli- 
sion models. Note though, that muons which land in 



an specific region of the material will sense fluctuating 
fields with a distribution characterized by a local A; and 
muons which land in other regions will sense fluctuating 
local fields with a different local A. Mathematically this 
amounts to making the time average (applying the strong 
collision model) before the spatial average (average over 
A with Eq. |4j 6 - 9 . /iSR signals with this Fluctuating 
Lorentzian Kubo-Toyabe (FLKT) model are shown in 
the upper panel in Fig. [T]for different fluctuation rates 
(vf). We want to notice two things. One is the strong 
depolarizing effect of dynamic fields (the polarization of 
the muon ensemble goes to zero at long times, while it 
is equal to 1/3 in the static case); and second, the de- 
crease in the relaxation rate of the polarization signal at 
high fluctuation rates. This last effect is known as the 
narrowing effect, and it happens because the local field 
is fluctuating so fast that the muons do not have time to 
precess much in the instantaneous local field. An impor- 
tant characteristic of the signal in the narrowing limit is 
its square-root exponential trend^^: 



A(i) = cxp(-4a 2 i/^/ 



,1/2 



(7) 



Now, we will assume that the magnetic moments are 
static but muons can hop from site to site. To account 
for this dynamical effect we have used the strong colli- 
sion model again , but now we have performed the time 
average after the spatial average. This means that a sin- 
gle muon can probe the whole distribution in Eq. [5] as it 
hops through the lattice (assuming again that the field 
from site to site is uncorrelated), and not only one spe- 
cific A. Then, muons are able to sense a bigger range of 
fields as they diffuse faster and faster. The lower panel in 
Fig. [1] shows several muon polarization signals using this 
Hopping Lorentzian Kubo Toyabe (HLKT) model. No- 
tice that no narrowing effect is observed. Further more, 
at high hopping rates the signal is exponential and inde- 
pendent of the hopping rat o 10 ! 11 : 



A(t) = exp(-4crf/3) 



(8) 



The absence of the narrowing effect in the HLKT 
model is a problem of the Lorentzian field distribution, 
which allows for the muon to sense unphysically large 
magnetic fields^. To appreciate this, we have calculated 
the probability that a muon senses a field bigger than 
Vh/l when hoping into a new site (at these fields the 
muons will be able to make a big precession before hop- 
ing to an other site and therefor reduce the polarization 
of the ensemble). At high hoping rates (vh 3> a) this 
probability is Aa/irvh- This can be used to calculate the 
probability that within the signal time window (typically 
10 /lis) the muon sense a highly depolarizing field. This 
probability is equal to ATa/it (T is the signal time win- 
dow) and it is independent of Vh- This means that as Vh 
is increased more and more the signal will tend to a Vh 
independent relaxation rate. This unphysical behavior 
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FIG. 1: Simulated signals for the FLKT model with a=l//s~ 
and at different fluctuation rates of the local magnetic field 
(in ^is _1 , upper panel); and for the HLKT model at different 
muon hoping rates (lower panel). The FLKT model shows 
the narrowing effect for IO^s" 1 and 100/xs" 1 , while the HLKT 
does not. 



of the HLKT model could be avoided by introducing a 
maximum field cutoff in the distribution on Eq. [SJ This 
cutoff is system dependent and should be calculated care- 
fully. This type of refinement of the model is outside the 
scope of this document. Never the less, the HLKT model 
is valid as long as v/, < a. 

We have extended the Lorentzian Kubo-Toyabe for- 
malism to account for systems where the magnetic mo- 
ments fluctuate and also muons diffuse. In this model, 
which we will call the Extended Lorentzian Kubo-Toyabe 
(ELKT) model, the dynamical behavior is also intro- 
duced by the strong collision model, but this time before 
and after the spatial average to account for the magnetic 
moment fluctuations and muon hoping respectively. We 
have studied this model numerically using a Monte Carlo 
algorithm. At time zero, this algorithm selects a A for 
the local field distribution with a probability given by 
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FIG. 2: Simulated signals for the ELKT model with a=l/is 
and Uf = 50^is and for different muon hopping rates (in /xs). 
The relaxation of the signal increases continuously with 
and the signal changes from a square root exponential to a 
simple exponential continuously. 



Eq. 21 Then, the local field is allowed to to fluctuate 
with a fluctuation rate Vf. The local field fluctuations 
are carried on until a muon-hop event happens, and this 
occurs with a frequency given by the muon hopping rate 
Vh- At this point a new A is randomly chosen from the 
probability distribution on Eq. |4] and, again, the local 
field is allowed to fluctuate with a field distribution char- 
acterized by the new field width. In between all these 
"collisions" the muon polarization function is properly 
evolved, and it is recorded into an array at a discrete set 
of times. 

Since our interest is to study the effect of muon diffu- 
sion in fast fluctuating paramagnetic materials (see Sec- 
tion HVJ, we have performed all our simulations in the 
narrowing limit (ff 3> a). Fig. [2] shows the effect of in- 
creasing the muon hopping rate on the /zSR signal from 
a diluted paramagnetic system. The signal is a root ex- 
ponential in the absence of muon diffusion and equal to 
Eq. [7] As the muon hopping rate is increased, the sig- 
nal shape changes continuously to an exponential at high 
diffusion rates. In the high hopping rate limit the signal 
is given by Eq. [5] 



III. PHENOMENOLOGICAL DESCRIPTION OF 
THE ELKT MODEL 

In order to use the numerical solutions of the ELKT 
model to analyze experimental data (see next section), we 
have fitted our simulated signals in the narrowing limit 
(vf S> a) to a stretched exponential function: 
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FIG. 3: Values of /3 (squares, left y-axis) and A (circles, right 
y-axis) as a function of vu- These values were found from 
fitting the simulated data, with a — 1/xs -1 and vj = 50/xs -1 , 
to Eq. [9] Lines are the plots of Eq. [10] (continuous line), [TT1 
(dotted line) and 1 121 (dashed line). 



A(t) = exp(-(\tf) 



(9) 



Fig. [3] shows the dependence of (3 and A on the muon 
hopping rate. Note that, as the muon hopping rate is 
increased from zero, the parameter (3 grows from 0.5 to 
1 as expected. We have found that the dependence of (3 
and A on the field fluctuating rate, muon hopping rate 
and a is well described by (see Fig. [3]): 
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(11) 



(12) 



where \i ow and Xhigh refer to different ranges where each 
equation apply. A; oto describes the Vh dependence for 
v h < 8a, and X h i g h for v h > 8a. 

We want to mention that, as a log(log(^4)) vs. log(t) 
plot shows, the numerical solutions of the ELKT model 
are close to stretched exponentials but are not equal. 
Nevertheless we have used this function because is com- 
monly used in the analysis of /iSR signals and it will 



allows us to understand previous published works (see 
Section [V}. In the next section we use Eq. [TU1 and [TTl to 
obtain quantitative measurements of the muon diffusion 
rate in Dyo.4Yi. gT^C^. 



IV. MUON DIFFUSION IN Dyo.4Y1.eTi2 Or 

The low temperature behavior of the ionic insulator 
Dy2Ti2C>7 have had much attention from researchers^. 
Below 20K, the Dy 3+ magnetic ions have an Ising charac- 
ter and they form a magnetically frustrated system due 
to their effectively antiferromagnetic interactions and to 
their arrangement in a lattice of corner shearing tetra- 
hedra. Below ~4K, the magnetic ions arrange such that 
two spins point inside every tetrahedra and two out. This 
arrangement of spins can be mapped into that of hydro- 
gen atoms in watter ice, and therefore this ground state is 
often referred as spin ice. More recently, it was noticed 
that the spin ice state is a background for topological 
excitations which behave as magnetic monopole a 13 ' 14 . 

We have used ^SR to study the effect of magnetic 
dilution in this system by substituting magnetic Dy 3+ 
ions with non magnetic Y 3+ . A single crystal of 
Dyo.4Yi.eTi207 was grown using an optical floating zone 
furnace. The single crystal was aligned using a Laue 
camera and cut such that the muons are implanted with 
their spin parallel to the (110) direction. All the sample 
preparation was performed in the Laboratory for Devel- 
opments and Methods in the Paul Scherrer Institut (PSI, 
Switzerland) . The /iSR measurements were performed at 
the GPS instrument at PSI, and which is located in the 
7rM3 beam line. We performed temperature scans be- 
tween 2K and 250K in an applied field of 100G transverse 
to the muon polarization. Also, we have performed com- 
plementary measurements with the field applied along 
the muon spin direction and in zero field. 

Fig. S] shows a ^SR signal from Dy .4Yi. 6 Ti 2 O7 at 
160K where the system is paramagnetic. Note in this 
picture the stretched exponential-like envelope of the sig- 
nal. For this reason, we first fitted the data to a cosine 
function times a stretched exponential envelope. The left 
panel in Fig. [5]shows that the parameter (3 increased with 
temperature from a value of «0.59 below 100K to «0.94 
at 250K, and A decreased monotonically from 21/is~ 1 to 
0.2/is -1 in the same range. As we mentioned in the pre- 
vious section, the ELKT model can explain the trend in 
f3 in terms of an onset of muon diffusion upon warming. 
Nevertheless, and according to Fig. [3J an onset of muon 
diffusion would produce an increase of A with tempera- 
ture. This is in contrast with what we observe. As we will 
see later, the decrease of A with temperature is also cap- 
tured by the ELKT model, and it is due to the fact that 
the thermally activated fluctuations of the Dy 3 + mag- 
netic moments are much bigger than the muon hopping 
rate (see Eq. [TT]) . Then, we analyzed our data by fitting 
the relaxing envelop to a stretched exponential where the 
relaxation rate is given by Eq. [TT]and the coefficient by a 
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FIG. 4: ^SR signal from Dyo.4Yi. 6 Ti 2 C>7 at 160K and in a 
transverse field of 100G. The continuous line is a fit to Eq. [9] 
multiplied by the experimental asymmetry (amplitude) and a 
cosine term. The dashed lines are the fitted envelope of the 
signal (i.e. without the cosine term). 

modified version of Eq. [TO] This modification consist on 
substituting the 1/2 factor on the right hand side of Eq. 
ITUlbv 0.41, such that /3 for T<120K is equal to the ex- 
perimentally observed value of 0.59 . This apparently ad 
hoc modification is justified by the fact that, even if the 
dilution is significant in our system, it is not in the very 
high dilution regime which the ELKT model assumes 
(20% of the DY/Y sites are occupied by magnetic ions in 
Dyo.4Yi. 6 Ti 2 07). As shown by Noakes-, as a system is 
magnetically diluted the relaxation function changes con- 
tinuously from a Gaussian Kubo-Toyabe to a Lorentzian 
Kubo-Toyabe. This means that in the narrowing limit, 
the coefficient of the stretched exponential will change 
from 1 to 0.5 continuously as the system is magnetically 
diluted. Our system is concentrated enough, such that 
this coefficient is slightly above 0.5 

It is well known that in Dy 2^207 the characteristic 
fluctuation rate of the magnetic moments in the param- 
agnetic state is thermally activate d 15 i 16 , and therefor the 
fluctuation rate of the local field at the muon site will be 
given by 6 : 

v f = v exp(-E/k B T) (13) 

where E is the energy gap between the ground state dou- 
blet and the first excited crystal field level. The reported 
values for this number ranges from 200K to 380 K 15 i 17 . In 
the absence of spectroscopic data, we used a value equal 
to 300K. In order to determine Uq we used a longitudinal 
field measurement at 60K to determine the fluctuation 
rate of the local magnetic field at this temperature. As- 
suming that the muons do not diffuse at these tempera- 
ture (this is indicated by a temperature independent (3 of 
approximately 0.59), we have obtained a fluctuation rate 
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FIG. 5: Left panel: coefficient of the power exponential re- 
laxing envelope (/3) and relaxation rate of the signal (A), from 
fitting the data of Dyo.4Yi.6Ti2 07 in zero external field with 
Eq. [9] times a cosine function. Right panel: muon hopping 
rate (vh) and characteristic width of the field distribution (a) 
as a function of temperature; obtained by fitting the data 
with A and /3 given by Eq. 1101 and [TT1 respectively. 

of 990MHz at this temperature. Using this value in Eq. 
ITS! we have obtained z^o=150GrIz, which is close to that 
found by neutron spin echo measurements (220GHz) 16 . 
Then, when we analyzed our data we set Vf to the value 
given by Eq. [13] and extracted a and Vh from the fit at 
each temperature. 

The fitted values for Vh and a are shown in the right 
panel in Fig. [5] Note the big error bars on Vh for T<80K. 
This only reflects the very steep grow of j3 at low values 
of Vh (see Eq. [TO)) . Figure [5] shows that muons start to 
diffuse between 120K and 140K, and the diffusion rate 
grows until «200K and then stays constant. Below 80K, 
the parameter a is constant and equal to lOO/is -1 . This 
is close to llO/is -1 , value that we found at 2K where 
the system is almost static within our time window. As 
the temperature grows the value of a decreases due to 
the different orientations that the magnetic moments can 
have in the higher crystal field level. 

V. COMMENTS 

As shown in the previous section, the ELKT model 
is able to explain the increase of /? with temperature in 
Dyo.4Yi.6Ti2C>7 in terms of the onset of muon diffusion. 
This increase of /3 with temperature has been observed in 
an other ionic insulato r 18 i 19 , but no explanation of this ef- 
fect was provided. Campbell et al?° observed an increase 
of /3 with temperature in the canonical spin-glass ^4(?Mn 
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at high temperature^. The authors mentioned that this 
effect could be produced by muon diffusion but no fur- 
ther analysis of the data was performed. This study also 
pointed out that the increase in j3 was accompanied by 
an increase, and subsequent saturation, of the relaxation 
rate of the signal upon cooling. The authors interpreted 
this as a consequence of muon retraping near the mag- 
netic impurities. In our case the relaxation rate of the 
signal decreases monotonically with temperature up to 
250K. Nevertheless, the ELKT model shows a saturation 
of v h above 160K (see Fig. [5]). It is not possible to know 
with the data taken if this saturation is due to muon re- 
traping (the slightly different ionic radius of Y and Dy 
create local distortions in the lattice^ which might fa- 
vor some muon sites over others). Nevertheless, we want 
to note that if a muon retraping process is present in 
Dyo.4Yi.6Ti207, then the ELKT model would provide 
only with an "averaged" muon diffusion rate since this 
model posses only one diffusion parameter (instead of the 
three dynamical parameters needed to describe diffusion 
with a single retraping mechanism^). Still, and even if 
muon retraping is present, the increase in /? from 0.59 to 
0.94 remains a signature of the onset of muon diffusion 
since its change is produced by the muon exploring dif- 
ferent regions of the material with different magnetic mo- 
ment densities (i.e. different characteristic field widths 
A). This fact is confirmed by our low transverse field 
measurements in a sample of Dyi.gYo.4Ti207, where we 
found that the parameter /3 is temperature independent 
and equal to 1 within error bars. This is expected by 
the theory since this sample is magnetically dense, and 
therefore all muon sites have approximately same local 
density of magnetic moments (i.e. the same field width 
A). 

We want to finish this section by reviewing a couple 
of limitations of the ELKT model. One limitation is the 
spherical nature of the internal field distribution which 



the ELKT model assumes (see Eq. [T]). Care should be 
taken when using Kubo-Toyabe models to fit single crys- 
tal data since the internal field distributions have not, 
in general, spherical symmetry. Nevertheless the biggest 
discrepancies are expected to be observed at low values of 
Vf and not so much at the high fluctuation rates explored 
in this work 18 . An other limitation of the ELKT is that 
it assumes that the magnetic fluctuations are character- 
ized by a single fluctuation rate of the local magnetic field 
(i>f). Then, it should not be used to model magnetic sys- 
tems with a richer dynamic response such as spin-glasses 
relatively close to T g . 



VI. CONCLUSIONS 

We have introduced the ELKT model which is suited 
to analyze /iSR signals from diluted magnetic systems 
where muon diffusion is present and where the mag- 
netic moments are dynamic. Also, we have provided a 
set of phenomenological formulas to extract muon diffu- 
sion rates from systems where the magnetic background 
fluctuates fast. We have used these formulas to analyze 
data from Dy .4Yi. 6 Ti 2 O7, and found that muons dif- 
fuse for T>140K. Also, we have found that the diffu- 
sion rate saturates for T>180K. We have shown that the 
increase with temperature of the parameter /3 in a hot 
diluted paramagnet is produced by the onset of muon 
diffusion, and this is further confirmed by our observa- 
tion that this parameter is temperature independent in 
a magnetically dense system. We believe that the use of 
the ELKT model can allow to study muon diffusion by 
introducing a small amount of strong magnetic impuri- 
ties in a system where muon diffusion can not be studied 
otherwise (e.g. in systems with small nuclear magnetic 
moments). 
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